Master Equation for Quantum Brownian Motion Derived by Stochastic Methods 
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The master equation for a linear open quantum system in a general environment is derived using 
a stochastic approach. This is an alternative derivation to that of Hu, Paz and Zhang, which was 
I ' based on the direct computation of path integrals, or to that of Halliwell and Yu, based on the 

, evolution of the Wigner function for a linear closed quantum system. We first show by using the 

' influence functinal formalism that the reduced Wigner function for the open system coincides with 

^N) I a distribution function resulting from averaging both over the initial conditions and the stochastic 

X I ■ source of a formal Langevin equation. The master equation for the reduced Wigner function can 

' then be deduced as a Fokker-Planck equation obtained from the formal Langevin equation. 

Si 

CN ■ I. INTRODUCTION 

\ Open quantum systems are of interest in condensed matter physics , quantum optics |^ , quantum measurement 
^ ■ theory | , nonequilibrium field theory (|-|^ , quantum cosmology ||9|]l0|] and semiclassical gravity |Tl|JT^ ] . An open 
C ■ \ quantum system consists of a subset of degrees of freedom, whose dynamics one is interested in, within a larger 
closed quantum system undergoing unitary evolution . This subsystem of interest is simply called the "system" 
whereas the remaining degrees of freedom constitute the "environment" . In general, the evolution of the system will 
be nonunitary and even non-Mar kovian. 

A typical example of an open quantum system is the quantum Brownian motion (QBM) model, which consists of 
\ a single massive particle interacting with an infinite set of independent harmonic oscillators with a Gaussian initial 
' state [ p^ . The coupling may be linear both in the system and environment variables or may be nonlinear in some 
or all of these variables. The frequencies of the environment oscillators are distributed according to a prescribed 
I • spectral density function, the simplest case corresponding to the so-called ohmic environment. Part of the interest 
^jpj' of the linear systems is that they are in many cases exactly solvable and detailed studies of different aspects of 
• • . open quantum systems can be performed. One of the issues that have received much attention in recent years is 
^ ' environment-induced decoherence as a mechanism to understand the transition from the quantum to the classical 
V : regime 

5^ Certain useful information for an open quantum system is contained in the master equation. The master equation 

is a differential equation for the time evolution of the reduced density matrix of an open quantum system. The master 
equation for linear coupling and ohmic environment at high temperature was first deduced by Caldeira and Leggett 
pTt , it was extended to arbitrary temperature by Unruh and Zurek |Q, and it was finally obtained for a g eneral 
environment (i.e. for an arbitrary spectral density function) by Hu, Paz and Zhang using path integrals [ [l9[ . This 
result was then extended to the case of nonlinear coupling by treating the interaction perturbatively up to quadratic 
order 

The reduced Wigner function is defined from the reduced density matrix by an integral transform [ pT|j22) . This 
function is similar in many aspects to a distribution function in phase space, although it is not necessarily positive 
definite, and the dynamical equation it satisfies is similar to the Fokker-Planck equation for classical statistical systems 
psj-psf. This equation is, of course, entirely equivalent to the master equation for the reduced density matrix and we 
will often also refer to it as the master equation. Halliwell and Yu exploited the fact that the Wigner function for a 
linear closed quantum system evolves according to the classical equations of motion to obtain the equation satisfied 
by the reduced Wigner function P6| . The reduced density matrix has been used to study decoherence induced by the 
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environment |27-31 



The Wigner function has also been used in studies of emergence of classicahty induced 



by an environment [32|, especially in quantum cosmology [|9|,|lO[ 



Langevin type of equations |3^,^ as a suitable tool to study the semiclassical limit have been used recently 
in semiclassical gravity and cosmology |35- |3^ , [l^ . In inflationary cosmology they have been used to describe the 
stochastic effect on the inflaton field |3^- 4£ ] or the stochastic behavior of large-scale gravitational perturbations , 
which is important for cosmological structure formation. So far, in the functional approach the Langevin equation 
has been mainly restricted to describe the classical or semiclassical limit. See, however, ref. pTf for a quantum version 
of the Langevin equation in operator language. 

A closer look at the influence functional, nevertheless, reveals that a formal Langevin equation can be extracted 
from this functional independently of the existence of a classical limit at least for quadratic influence actions. This 
Langevin equation is used to show that the reduced Wigner function can be written as a formal phase-space distribution 
function associated to a stochastic process |Q (as earlier suggested in ref. |^^). The master equation governing its 
time evolution is then deduced as the corresponding Fokker-Planck equation. 

The plan of the paper is the following. In sec. O we briefly summarize the essential concepts and results of the 



influence functional formalism for linear open quantum systems. In sec. Ill we show how the reduced Wigner function 



for the swtem can be expressed as an average over the different realizations of a stochastic process. This result is used 



in sec. IV to give an alternative derivation of the master equation for a general environment. Finally, we summarize 
and discuss our results in sec. 0. 

II. INFLUENCE FUNCTIONAL FORMALISM AND MASTER EQUATION FOR LINEAR OPEN 

QUANTUM SYSTEMS 

Let us first review a QBM model as an example of linear open quantum system. We consider a harmonic oscillator 
of mass M, the "system" , coupled to a bath of independent harmonic oscillators of mass m, the "environment" . For 
simplicity, let us assume that the system and environment are linearly coupled. The action for the whole set of degrees 
of freedom is defined by: 

S[x, {q,}] = S[x] + S[{q,}] + S,^t[x, {q,}], (2.1) 

where the terms on the right-hand side correspond to the action of the system, the environment and the interaction 
term respectively. They are given by: 

S[x] = J dt{]^Mx^ - iikTf^^x^), (2.2) 
S[{1,]] ^Y.j d^^l"^^' - l^^h')^ (2-3) 



Sint[x,{qj}] ^y^Cj [ dtx{t)qj{t) = ( duj ^"^^ I{uj) ( dtx{t)q{t; 
, J Jo Trc(uj) J 



to), (2.4) 



where we introduced the spectral density I{uj) — i:c^j{2mijjj) ^5{ui — Uj) in the last equality, c{uj) and q{t;uj) are 
functions such that c{ujj) = Cj and q{t;ujj) = qj{t), Cj being system-environment coupling parameters, and O and 
LUj are, respectively, the system and environment oscillator frequencies. When no special form is assumed for the 
spectral density /(w), this is usually referred to as a general environment. One of the most common particular cases 
is the so-called Ohmic environment, defined by I{lo) ~ ui (some high frequency cut-off may be sometimes naturally 
introduced) . 

The reduced density matrix for an open quantum system is defined from the density matrix p of the whole system 
by tracing out the environment degrees of freedom 

Pr{xf,x'j:,tf)^ Y[dqjp{xf,{qj},x'f,{qj},tf) = dx^dx'^J{xf,x'j:,tf;x,,x'^,U)pr{x^,x'^,U), (2.5) 
J ^ J 

where the last equation gives the evolution of the reduced density matrix by means of the propagator J, which is 
defined in a path integral representation by 
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J{xf,x'f,tf-x,,x[,U)^ J Vx J Pa;'e^(^[^l-^[^'l+^^-[^•^'l)/^ (2.6) 

x(ti) — Xi x' {ti)—x'. 

where Sif[x, x'\ is the influence action introduced by Feynman and Vernon [^o| . When the system and the environment 
are initially uncorrelated, i.e., when the initial density matrix factorizes (piU) = pr{ti)(E> Pe{ti), where Pr{ii) and pdU) 
mean, respectively, the density matrix operators of the system and the environment at the initial time) the influence 
functional, defined by F[x,x'] — exp^iSipix, x'])/h, can be expressed in the following way: 



'jr{S[{q,}]^S[{q'^}] + S[x,{qj}] 



F[x,x']=l[l dq^dqfdqf J Vq, J Pq^ exp 

},t,). (2.7) 

When the initial density matrix for the environment pe{{qj^^}, {qf^''}, U) is Gaussian, the path integrals can be exactly 
performed and one obtains ||50| , |l7|] : 

Sif[x, x'] = -2 ^ ' ds ^ ds'A{s)D{s, s')X{s') + ^J^''^^J^' ds'A{s)N{s, s')A{s'), (2.8) 

where X{s) = {x{s) + x'{s))/2 and A(s) = x'{s) — x{s). The kernels D{s, s') and N{s, s') are called the dissipation 
and noise kernel, respectively. 

For environments consisting of a n infinite number of oscillators it is especially convenient to rewrite the first term 
on the right-hand side of Eq. ( |2.S| ) as 

' dsj^' ds'A{s)Htare{s, s')X{s'), (2.9) 

where we defined Hi,are{s, s') as formally equivalent to — 2Z?(s, s')9{.s — s'). Being the product of two distributions the 
latter expression is not well defined in general and suitable regularization and renormalization may be required; see 
pl[ for details. The local divergences present in Hbareis,s') = H{s,s') + Hdiv5{s — s') can be canceled by suitable 
counterterms fldiv in the bare frequency of the system O = Qren + ^div From now on we will consider that this infinite 
renormalizatio n, if necess ary , has already been performed so that both Qren and H{s, s') are free of divergences. 

From Eqs. (^^) and (|]^) a differential equation for the system's reduced density matrix known as the master 
equation can be derived. The expression for a general environment was first obtained by Hu, Paz and Zhang using a 
path integral approach (see for a slightly different derivation): 

-!fiA(i)(i - i') ('J^ - ^) ft, + hB(t){x - i') (^ + A'j _ iMC{t){x - x'fpr, (2.10) 

where the functions (5il^(t), A{t), B{t) and C{t) represent a frequency shift, a dissipation factor and two diffusive 
factors, respectively. For explicit expressions of these functions see sec. IV. An alternative representation for the 
system reduced density matrix is the reduced Wigner function Wr {X, p, t) defined as 

1 

Wr{X,p,t)^— dAe'P'^/^priX - A/2,X + A/2,t). (2.11) 
27rn 



It follows immediately that the master equation (2.10) can be written in the following equivalent form: 

^ - {H,, WApb + mtf-^ + ^5(t)^ + r.MC[tf^, (2.12) 

where {HR,Wr}pB = -{p/M)dWr/dq + Am\{t)qdWr / dp with n\{t) = Vi^ + 6VL^{t). This equation was directly 
derived by Halliwell and Yu exploiting the fact that the Wigner function for the whole closed quantum system 
evolves according to the classical equations of motion. Note that Eq. ( |2.12| ) is formally similar to the Fokker-Planck 
equation for a distribution function. 
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III. STOCHASTIC FORMAL EXPRESSION FOR THE REDUCED WIGNER FUNCTION 



In this section we show that the reduced Wigner function can be written as a formal distribution function for some 
stochastic process (see ref. for a detailed exposition). This will be the key starting point in the derivation of the 
master equation given in the next section. 

In order to find an explicit expression for the reduced density matrix ( |2.5| ) at a time tf, we need to compute the 
path integrals appearing in Eq. (^_^) for the reduced density matrix propagator. From now on we will consider h = 1. 
After integrating the system action by parts and performing the Gaussian path integral for A(t) with and Af 
fixed, we obtain 



where L{t, t') = M ^^tt + ^renj S(t — t') + H{t, t'). Taking into account the surface terms arising from the integration 
by parts of the system action and definition ( p.ll ) for the reduced Wigner function, the result of the integration over 



Ai gives 



PriXf-Af/2,Xf + Af/2,tf)=27ridct—] J dX, j ^ I?Xe"*(^-^)-^"(^-^)e-'*^^/^/iy,(X„ MX,, t,)- 

(3.2) 

The next step to perform is the following functional change: 

X{t) {x, = X{U), = MX, = MX{U), e(t) = [L ■ X)it) } . (3.3) 

Note that with this change the function X{t) gets substituted by the initial conditions {Xi,pi) and the function ^{t) 
in the path integration. It is important to note that at this point the function ^(i) is not a stochastic process but just 
a function over which a path integral is performed. The functional change (^^) is invertible as can be explicitly seen: 

{X„ p„ C(<)} X{t) = Xo{t)+ dt'Gret{t,t')^{t'), (3.4) 

where Gret{t',t") is the retarded (i.e., Gret{t',t") — for t' < t") Green function for the linear integro-differential 
operator associated to the kernel L{t,t'), and Xinhit) — J^.. dt'Gret{t,t')^{t') is a solution of the inhomogeneous 
equation {L-Xinh){t) — £,{t) with initial conditions Xinh{ti) = and dXinh{t') / dt'\^,^^, — 0. On the other hand, Xo{t) 

is a solution of the homogeneous equation {L-Xo){t) — 0, with initial conditions Xo{ti) = Xi and X{ti) — Pi/M . Since 
the change is linear, the Jacobian functional determinant will be a constant (this can be clearly seen by skeletonizing 
the path integral). After performing the functional change, we obtain 

/OO f'OO (> 

dX, / dp, / I?e'5(X(i/)-X/)e-^«-^-'-«e-'*^^(*/)'^/W^,(X„p„t,), 
OO J — OO J 

(3.5) 

where the delta function 5{X{tf) — Xf) was introduced to restrict the functional integral J with free ends, in order 
t o ta ke into account the restriction on the final points of the allowed paths for the integral /'^^ VX appearing in Eq. 



(3.2). The contribution from the Jacobian has been included in the constant K. By demanding the reduced density 
matrix to be normalized, i.e., that Trpr{tf) = 1, provided that the initial Wigner function is properly normalized, 
this constant can be determined to be 



K 



[det(27rX)] ^. (3.6) 



Finally, using the definition (2.11) for the Wigner function and the fact that {2tt) ^ J^^dAfC^'^f'^fe' 



-iMX{tf)A 



S{MX{tf) — pf), we get an expression for the reduced Wigner function 
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Wr{Xf,pf,tf)^K dxJ dpj V^S{X{tf)~Xf)S{MX{tf)~pf)e-^^-''~'<WriX,,p,,U), (3.7) 

J —oo J —oc J 

which can be written in the foUowing suggestive way: 

Wr{Xf,pf,tf) = (^(^5{X{tf)-Xf)S{MX{tf) -Pf))^)^ ' (3-8) 

where 

(...). = [det(2^iV)]"* / Pe-e"*^'^"''^ (3.9) 



dX, / dp,...Wr{X„p,,t,). (3.10) 

-oo J —oc 

Thus the reduced Wigner function can be interpreted as an average over a Gaussian stochastic process ^{t) with 
{^{t))^ ~ and {^{t)£,{t')) ^ — N{t,t') as well as an average over the initial conditions characterized by a distribution 
function Wr{Xi,pi,ti). It is only after formally interpreting ^{t) as a stochastic process characterized by Eq. ( |3.9| ) 
that the equation defining ^(t) in the functional change (|3.3[) 



{L-X){t)^at), (3.11) 



can be regarded as a Langevin equation. We insist that, in general, Eq. ( |3.11 ) is not meant to describe the actual 



trajectories of the system, but it should rather be regarded as a formal tool. We should also remark that X{tf) and 



X{tf) in Eq. ( |3.8[ ) correspond to solutions of the Langevin equation (3.11) for a given realization of ^(t), and that 
Xf and pf are coordinates of a point in phase space. 

Note, in addition, that although Wr{Xi,pi,ti) is real, which follows from the hermiticity of the density matrix, 
and properly normalized, in general it is not positive everywhere and, thus, cannot be considered as a probability 
distribution. The fact that the Wigner function cannot be interpreted as a phase space probability density is crucial 
since most of the nonclassical features of the quantum state are tightly related to the Wigner function having negative 
values. For instance, a coherent superposition state is typically characterized by the Wigner function presenting strong 
oscillations with negative values in the minima ||3^ , |3l|] , which are closely connected to interference terms. 

IV. FROM LANGEVIN TO FOKKER-PLANCK: DERIVATION OF THE MASTER EQUATION 

As mentioned above there is a sim ple one-to-one correspondence between any density matrix and the associated 



Wigner function introduced in (2.11). Taking this correspondence into account, the equation satisfied by the reduced 
Wigner function is equivalent to the master equation satisfied by the reduced density matrix. 

Equation (pM shows that the reduced Wigner function can be interpreted as a f ormal distribution in phase space. 



By deriving it with respect to time and using the Langevin- type equation in ( |3.11 ), one can obtain a Fokker-Planck 



differential equation describing the time evolution of the system's reduced Wigner function. 

The derivation of the Fokker-Planck equation from the Langevin equation with local dissipation is well understood 
(see ref. p^). However, in our case the existence of n onlo cal dissipation makes it convenient to review the main steps. 



Let us begin by computing dWr/dt from expression (3.8) 
dWriX,p,t) 



dt 



= Ux{t)5'(X{t)-X)5{MX{t)-p)) \ + IU{X{t)~X)MX(t)6'{MX{t)-p)\ \ 

\ ^/ Xi,pi \ ^/ Xi.pi 

" ' ({siXit)-X)MXmMXit)-p)) ) , (4.1) 



M dX dp 



where the fact that X{t), d/dXit) and d/dX{t) may be replaced by p/M, —d/dX and —d/dp respectively, since they 
are multiplying the delta functions, was used in the second equality. L et u s now concentrate on the expectation value 
appearing in the last term and recall the expectation values defined in ( |3.9[ )-( pTTo| ). We will consider the Langevin- type 
equation 

{L-X)it')=m, (4.2) 
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corresponding to the func tional change ( |3.3[ ) and substitute the corresponding expression for MX(t) so that the last 
expectation value in (O) can be written as 



rt 

ren 



Mnl^„XWr{X,p,t) + (((-/ dtH{t,t')X{t')+at) ] 6{X{t)-X)5{MX{t) -p)) ) . (4.3) 



Any solution of Eq. (|4.2|) can be written as 



X{t') - X,(t') + / dt"Gad.{t'XW), (4.4) 



where Xh{t') is a solution of the homogeneous equation {L ■ X){t') — such that Xh{t) = X, Xh{t) = p/M and 
Gadv{t' 1 1") is the advanced {i.e., Gadv{t' , t") = for t' > t") Green function for the li nea r integro-differential operator 
associated to the kernel L{t,t'). The particular solution of the inhomogeneous Eq. (|4.2|) 



X,nh{t') = dt"Gadv{t',t")at") (4.5) 

has boundary conditions Xinh{t) — 0, dXinh{t')/dt' — 0. Both Xfi{t') and Gadv{t' ,t") can be expressed in 

t'=t 

terms of the homogeneous solutions ui(<') and U2{t'), which satisfy ui{ti) = 1, ui[t) = and ^2(^1) = 0, U2(i) = 1 
respectively: 

Xu{t') ^ X U{t') ^-i(t')) + ^-i(^'), (4.6) 
V ui[t) J Ul(t) 

r u: 1 U^{t')u2{t")-U2{t')u^{t") 

Gadv[t ,t ) ^ - — — — — 7-JJTtl(t -t). (4./'j 

M Ui(t")U2[t") — U2{t")Ui[t") 

We use the advanced propagator so that there is no dependence on the initial conditions at time t' — U coming from 
the homogeneous solution but just on the final conditions at time t' = t, i.e., on those the Fokker-Planck equ ation is 
written in terms of. Using expression (4^) the first term within the expectation value appearing in Eq. (^^) can be 
reexpressed as 

dtH{t, t') ( (x{t')S{X{t) - X)5{MX{t) - p) 

dt'H{t,t')Xhit')Wr{X,p,t) + dt' j^^ dt"Hit,t')Gadvit',t")(^(^at")HX{t)-X)SiMXit)-p))^^ . (4.8) 

The first term on the right-hand side can in turn be written as 

- {M5n{t)X + 2A{t)p) Wr{X,p,t), (4.9) 

where 

1 r* 

sm = — dt'H{t,t')[u2{t') - {u2{t)/Mt))Mt')], (4.10) 

A(i) = i(Mui(i))-i f dt'H(t,t')u^{t'). (4.11) 

In order to find an expression for {^^{t')5{X{t) — X)5{MX(t) — p)^ we use Novikov's formula for Gaussian stochas- 
tic processes |5^, which corresponds essentially to use (|3.9| ) and functionally integrate by parts with respect to £,{t), 

{mm e])^ = /* dt"N{t', t") {SF{t- ^]/san)e ■ (4.12) 

We then obtain the following expression: 



ti 
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dt" 



dt"N{t', t") 



+ 



5X{t"') 6 
d^{t") 5X{t"') 



\ S^{t") 6X{t"') 
S{X{t) - X)5{MX{t) -p) 



dt'" / dt"N{t', t")5{t"' -t)(- 



6{X{t) - X)S{MX{t) -p) 



5X{t"') d , j^j^^^i'") 9 



S^it") dX 6^{t") dp 



^ \5i{t")dX^ 5i{t")dp^ 



5{X{t) - X)6{MX{t) -p)(_ 



(4.13) 



where we used again the presence of the delta functions to substitute the functional derivatives S/SX{t"') and S/dX{t"') 
by —5{t"' — t) ■ d/dX and —5{t"' ~t) ■ M ■ d/dp, respectively, in the second equality. Functionally differentiating with 
respect to £,{t") expression (Oh for X{t) and analogously for X{t) we get 



5X{t') 

m") 



Gret{t' , t"), 



which after substitution into ( f.l3| ) leads to 



^{t')SiXit) - X)SiMXit) - p)) ^ = - / dt"N{t',t") ( Gret{tX)Qj^+M 



d 



dt' dp 



(4.14a) 
(4.14b) 

Wr{X,p,t). (4.15) 



The retarded Green function can also be expressed in terms of the solutions of the homogeneous equation ui{t) and 
U2{t), which were previously introduced, as 

r 1 ui{t')u2{t") - U2it')ui{t") 

Gret(t ,t ) ^ 17^7777^ — 7W\ w777^ — 7777T^(* ~ * ■ (4.16) 

M Ui{t")U2{t") - U2it")ui[t") 

Note that it is important to use now the expression in terms of the retarded propagator Gret and the initial conditions 
Xi and Pi (at time t' = ti), since the "fina l" conditions X(t) and MX{t) depe nd o n £,{t") (for t" < t). Putting all the 
terms together, i.e., (4.3), (4.8) and (4.15), we reach the final expression for (4.1): 



dt 



dp 



dXdp 



dp"^ 



where 6Q.{t) and A(t) are given by Eqs. (4.1C) and (4.11), and 



B{t)= / dt"'N{t,t"')Gret{t.t"')- / dt'H{t,t') / dt"GaUt',t") / dt'" N (t" ,t"')Gret{t.t"') 



C{t) = / dt"'N{t, t"') 



dt 



dt'H{t,t') / dt"Gadv{t',t") / dt'"N{t",t'") 



dGret {t, t" 



df' 



(4.17) 

(4.18) 
(4.19) 



The last two expressions were obtained by combinin g th e second term within the expectation value appearing in (4.3) 
and the second term on the right-hand side of Eq. ([4.81) . It sh ould be taken into account that if we put back the ?i's, 
there appears one with every noise kernel in Eqs. ( [4.1^ ) and (4.19). 

The expressions ( |4.10| ), ( [4.11D , ( p^ ) and ( [4.19D for Sn{t), A{t), B{t) and C{t), respectively, coincide exactly with 
those of ref. [^ , which are in turn equivalent to those obtained in ref. [p^ . Thus, this derivation of the master 
equation based on a stochastic description for the system is an alternative to those given previously |2^,|l^,^ and is, 
of course, in agreement with their results. 
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V. DISCUSSION 



In this paper we have considered the stochastic description of a hncar open quantum system to give an alternative 
derivation of the corresponding master equation. We have shown that the reduced Wigner function can be written as a 
formal distribution function for a stochastic process characterized by a Langevin-type equation. The master equation 
has then been deduced as the corresponding Fokker-Planck equation for the stochastic process. This derivation can 
be extended to the case of nonlinear interaction between system and environment by computing perturbatively the 
influence action up to quadratic order and even to the case of a general potential for the system [m] ] . 

It should be pointed out that whereas one can derive the Fokker-Planck equation from the Langevin equation, the 
opposite is not possible in general. One can always consider Langevin equations with stochastic sources characterized 
by different noise kernels which, nevertheless, lead to the same Fokker-Planck equation and, thus, the same master 
equation. This can be argued from the expressions obtained in the derivation of the Fokker-Planck equation. Let us 
consider, for simplicity, the situation corresponding to local dissipation. A local contribution to the noise gives no 
contribution to B{t), but it does contribute to C(t) as can be seen from Eqs. (4.18) and (4.19) taking into account 
that Gret{t,t) = and dGret{t' ,t)/dt'\^,^^ = M~^. Thus, one can always choose any noise kernel that gives the 
desired B{t) and then add the appropriate local contribution to the noise kernel to get the desired C{t) keeping B{t) 
fixed. Note that changing the noise kernel does not change A{t). To illustrate the fact that there exist different noise 
kernels giving the same B{t), as was stated above, one may consider the particular case corresponding to the weak 
dissipation limit so that Gret{t,t') ~ {MVt)-^ s\an{t - t')0{t - t'). To see that a different N{t,t') giving the same 
B{t) as N{t,t') exists reduces then to show that there is at least one nontrivial function v{s,t) = N{t,t') — N{t,t') 
(with s = t — t') such that for any t ds sin(r2s)i^(s, t) = 0, which can be shown to be the case. 

The fact that different Langevin equations lead to the same master equation reflects that the former contains more 
information than the latter. To be more precise, what we showed was that a Langevin equation contains in general 
more information that the corresponding Fokker-Planck equation. To extend this assertion to the master equation, 
one should make sure that different Langevin equations leading to the same Fokker-Planck equation can be obtained 
from an influence functional. Indeed this can be shown to be the case provided that one considers general Gaussian 
initial states for the environment. The inequivalence between the Langevin equation and the master equation can be 
qualitatively understood in the following way. In the influence functional it is only the evolution of the environment 
degrees of freedom that is traced out. Of course, having integrated over all the possible quantum histories for the 
environment, no correlations in the environment can be obtained. Nevertheless, since the system is interacting with 
the environment, non-Mar kovian correlations for the system at different times may in general persist. On the other 
hand, when considering either the reduced density matrix or its propagator, also the system evolution, except for 
the flnal state, is integrated out. Consequently, information on non-Markovian time correlations for the system is 
no longer available. Thus, only when the system's reduced dynamics is Markovian, i.e. the influence functional 
is local in time, we expect that the Langevin equation and the master equation contain the same information. In 
particular, for a Gaussian stochastic source, as in our case, the Langevin equation contains the information about 
the system correlations at different times, which the Fokker-Planck equation cannot in general account for. Only in 
the case in which the dynamics generated by the Langevin equation is Markovian one can compute the correlation 
functions just from the solutions of t he Fokker-Planck equation or, equivalently, the master equation for the propagator 
J{x2, X2,t2] xi, x'ljti); see Eq. (2.6). The key point is the fact that the propagator for the reduced density matrix 
only factorizes when the influence action is local. See ref. [^8|| for a detailed argument on this point. 

It is important to note that for a closed quantum system the evolution determined by the time evolution operators 
U{t2,ti) obtained from the Schrodinger equation is always unitary and, thus, also Markovian. That is why the 
Schrodinger equation suffices to get the correlation functions for a closed quantum system. On the contrary, for an 
open quantum system the evolution is nonunitary and, provided the influence action is nonlocal, not even Markovian. 

Finally, we should insist on the fact that, although we have exploited the formal description of open quantum 
systems in terms of stochastic processes, a classical statistical interpretation is not always possible. Thus, although 
the Wigner function is a real and properly normalized function providing a distribution for the initial conditions of our 
formal stochastic processes, it is not a true probability distribution function in the sense that it is not positive definite 
in general. In fact, this property is crucial for the existence of quantum coherence for the system. Nevertheless, 
even though the Langevin equation does not in general describe actual classical trajectories of the system, it is still a 
very useful tool to compute quantum correlation functions ||4^ ] or even as an intermediate step to derive the master 
equation. 
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